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Understanding the behavior of quantum systems subject to magnetic helds is of fundamental im¬ 
portance and underpins quantum technologies. However, modeling these systems is a complex task, 
because of many-body interactions and because many-body approaches such as density functional 
theory get complicated by the presence of a vector potential into the system Hamiltonian. We use 
the metric space approach to quantum mechanics to study the effects of varying the magnetic vector 
potential on quantum systems. The application of this technique to model systems in the ground 
state provides insight into the fundamental mapping at the core of current density functional theory, 
which relates the many-body wavefunction, particle density and paramagnetic current density. We 
show that the role of the paramagnetic current density in this relationship becomes crucial when 
considering states with different magnetic quantum numbers, m. Additionally, varying the magnetic 
held uncovers a richer complexity for the “band structure” present in ground state metric spaces, 
as compared to previous studies varying scalar potentials. The robust nature of the metric space 
approach is strengthened by demonstrating the gauge invariance of the related metric for the para¬ 
magnetic current density. We go beyond ground state properties and apply this approach to excited 
states. The results suggest that, under specihc conditions, a universal behavior may exist for the 
relationships between the physical quantities dehning the system. 
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I. INTRODUCTION 

Systems immersed in magnetic fields are a fundamen¬ 
tal research topic, as is the case, for example, for atoms 
immersed in strong fields mn], and are also an integral 
part of emerging quantum technologies, such as quantum 
computation, which utilise quantum systems controlled 
or otherwise affected by magnetic fields. For example 
the inhomogeneous magnetic field generated by the nu¬ 
clei’s spins decreases quantum coherence of electron spin 
qubits in III-V quantum dots [3], while full polariza¬ 
tion of the spin bath through an applied magnetic field 
suppresses electron-spin decoherence in nitrogen-vacancy 
centers and nitrogen impurities in diamond 0. Under¬ 
standing systems immersed in a magnetic field at a quan¬ 
tum level is therefore of both fundamental and techno¬ 
logical importance. 

However, the presence of a magnetic field introduces 
additional complexity to the system’s Hamiltonian. In 
fact, as opposed to a confining potential, which is de¬ 
fined by a scalar potential, V (r), the magnetic field is 
defined by B (r) = V x A(r), where A (r) is a vector 
potential. To account for its presence, density functional 
theory (DFT) must be extended to current density func¬ 
tional theory (CDFT). In this paper we will use the met¬ 
ric space approach to quantum mechanics [SHE] to study 
the effect on quantum systems of changing the vector 
potential and we will carefully consider the implications 
of the results for CDFT. This is particularly relevant as 
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there are still open questions with respect to the funda¬ 
mentals of this theory (see, for example. Refs. [5lH3j). 

The metric space approach involves the derivation of 
“natural” metrics from conservation laws to assign a dis¬ 
tance between two physically relevant functions [8]. In 
recent work these metrics were applied to the basic vari¬ 
ables of both standard DFT [MZllIl] and CDFT 0. In 
Refs. [ 5 ] and 0 it was demonstrated that the core the¬ 
orems of DFT and CDFT, respectively, indeed represent 
mappings between metric spaces: This helps understand¬ 
ing the power of the metric space approach and why its 
use has already allowed for the discovery of additional 
properties of these core theorems. The results in Ref. 0 , 
pertaining to systems with fixed magnetic fields, consid¬ 
ered only the effects of varying the scalar potential; in this 
paper we demonstrate how the metric space approach to 
quantum mechanics can be applied to analyse systems 
while varying the vector potential and hence the mag¬ 
netic field. 

A significant issue for any theory involving magnetic 
fields is gauge transformations of the scalar and vector 
potentials. The magnetic field, along with all physical ob¬ 
servables, is gauge invariant. Hence, in order to properly 
describe the distances between physical quantities, the 
metrics we derive must be robust against gauge transfor¬ 
mations. However, quantities such as the wavefunction 
and the paramagnetic current density are gauge variant 
and changes in the vector potential can result in gauge 
transformations for these quantities. Thus, here we will 
extend the metric space approach to ensure that the met¬ 
rics associated with these quantities are gauge invariant. 

We will provide further insight into the fundamental 
mappings between key physical quantities at the core of 
CDFT by studying the ground state of model systems as 
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the vector potential is varied. In particular we will ex¬ 
amine how the “band structure” introduced into ground- 
state metric spaces by the presence of a magnetic field [5] 
responds to changes in the field. To complement this pic¬ 
ture, we will apply the metric space approach to quantum 
mechanics to explore the properties of excited states. Re¬ 
sults will also help with validating the conclusions from 
the ground-state analysis. 

The rest of this paper is organised as follows. In Sec.[TT| 
we briefly review how functions obeying integral conser¬ 
vation laws can be cast as metric spaces and the ap¬ 
plication of this approach to the wavefunction, particle 


density, and paramagnetic current density. Section III 


demonstrates how the gauge properties of the paramag¬ 
netic current density are accounted for when forming the 
related metric space in order to ensure that the metric is 
gauge invariant. In Sec. IV we examine how the “band 
structure” present in metric spaces for ground states is 
affected by variations in the magnetic field and the rele¬ 
vance of this for CDFT. Section |V] goes beyond ground 
states. In Sec. VI we present a summary and our con¬ 
clusions. We use atomic units h = e = nif, = l/47reo = I 
throughout this paper. 


II. METRIC SPACES FOR PHYSICAL 
FUNCTIONS 

In Ref. [5] a general procedure for deriving metric 
spaces from conservation laws was presented. For com¬ 
pleteness in this section we briefly review this procedure 
and the properties of these metric spaces. 

A metric space {X, D) consists of a non-empty set of 
points X and a metric H : A x AT —>■ K, which assigns a 
distance between any two elements of X. For all a,b,c G 
X the metric must satisfy the following axioms nulls]: 


D {a, &) > 0 and D (a, b) = 0 a = b, (I) 
D{a,b) = D{b,a), (2) 

D (a, b) ^ D (a, c) + D (c, b). (3) 

These axioms are known as positivity, symmetry, and the 
triangle inequality respectively. 

Consider now a conservation law of the form 

J l/(x)rdx = c, (4) 

with c a finite, positive constant and x = {xi,... ,xn) 


including any spatial or spin coordinate in any dimen¬ 
sionality. Then for each p such that 1 ^ p < oo, the 
entire set of functions that satisfies Eq. Q forms an 
vector space. Then the corresponding metric 


Df (/l,/2) = 


l/l (x) - /2 (x)l^dx 


(5) 


also applies to the restricted set of physical functions 
obeying the conservation law @ 0. As c spans the set 



FIG. 1. (Color online) Sketch of the structure of the “onion- 

shell” geometry, consisting of a series of concentric spheres. 
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The first three spheres are shown, with radii c? ,i = 1, 2, 3. 


of its physically allowed values {ci}, the metric ([^ im¬ 
poses on its metric space an “onion-shell” geometry that 

I 

consists of a series of concentric spheres with radii cf , as 
sketched in Fig. 

We note that the procedure developed in Ref. [H] can 
be extended to conservation laws of the form 


= ( 6 ) 

i=l 


as the P vector spaces for sums are directly analogous to 
the LP spaces for integrals m- In this case the induced 
metric will be 


Df (/l,/2) 


n p 


j=i 


(7) 


Thus, we have a general procedure to construct a met¬ 
ric for any conservation law that is, or can be cast, in 
the form of Eq. Q or ([^. We can therefore state that 
such conservation laws induce metrics on the set of re¬ 
lated physical functions. As they descend directly from 
conservation laws, these “natural” metrics are non-trivial 
and contain the relevant physics. 
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Following this procedure the metrics 


= j + IV'2|^) dri... drN 

-2 J tpl^p2dri... dvN 

Dp{pi,p2) =y |pi(r)-p2(r)|dr, 

(jpi>jp 2 ) = J Kr X [jpi (r) - jp, (r)]}J dr, (10) 


( 8 ) 

(9) 


for wavefunctions '0 [5] , particle densities p [5] , and para¬ 
magnetic current densities jp = {jp^,jpy,jp^) [H] have 
been introduced. These metrics follow from, respectively, 
the conservation of wavefunction norm, of particle num¬ 
ber, and of the z component of the angular momentum; 
the latter can be expressed as Eq. when using the 
relation [5] 


y [rxjp(r)]^ dr = (0111^10) =771. (11) 

For paramagnetic current densities, Eq. ( |11[ ) directly im¬ 
poses an equivalence relation on the set of all param¬ 
agnetic current densities because [r x jp (r)]^ is indepen¬ 
dent of jp^. As a result, is a metric defined on a 
set of equivalence classes of paramagnetic current den¬ 
sities, with the classes characterised by paramagnetic 
current densities with the same transverse component 
jp7 ~ (jpa: ) jpv ) ■ 

The radii of the concentric spheres in the “onion-shell” 
geometry of the aforementioned metric spaces are: \/N 
for the wavefunction metric space, N for the parti¬ 
cle density metric space, and \m\ for the paramagnetic 
current density metric space When considering in¬ 

dividual spheres, the diameter of the sphere imposes an 
upper bound on the value of the distance. 


III. GAUGE INVARIANCE OF METRICS 

When dealing with electromagnetic fields, it is impor¬ 
tant to consider the choice of gauge. The scalar and 
vector potentials in the Hamiltonian are not unique, as a 
change of gauge transforms the potentials according to 

W (r) = U (r) + c', A'(r) = A (r) - Vx, (12) 

where c' is a constant and x (r) is a scalar field [18]. These 
transformations preserve the electromagnetic fields and 
all physical observables. 

With regard to the quantities we consider in this paper, 
the particle density is gauge invariant, but the wavefunc¬ 
tion and paramagnetic current density are not. After a 
change of gauge, the wavefunction undergoes a unitary 
transformation, which is given by m 

0'(r) = (r). (13) 


The paramagnetic current density transforms according 
to [iH] 

jp(r) =jp(r) + P(r)Vx- (14) 

Thus, when considering changes in the vector potential, 
we must be aware of the effect of gauge transformations 
on the physical quantities we are considering. Our met¬ 
rics are constructed to provide non-trivial information 
that is physically relevant; they are based in fact on con¬ 
servation laws. It is paramount then that they are also 
gauge invariant. The issue of ensuring that the wave- 
function metric is gauge invariant has been discussed in 
Ref. [Sj; we provide a formal review of the approach in 
Appendix]^ In this paper we wish to discuss gauge in¬ 
variance with respect to the paramagnetic current den¬ 
sity metric. 


A. Gauge invariance for the paramagnetic current 
density metric 


To consider the gauge properties of (jpi Jpa); 
of all we require that (r) and (r) are within the 
same gauge. Then, applying the gauge transforma¬ 
tion ([m]), we obtain 


D 


(jpi’jp2) [jpi W-jp2 W]}. 

= y 1 ( 1 - X {jp^ (r) - jp, (r) 


dr, 


+ [Pi (r) - P 2 (r)] VxD J dr. (15) 


Equation (15) states that, in general, the paramagnetic 


current density distance defined by Eq. (10) is modified 


by a gauge transformation. This seems to contradict the 
fact that we base our metrics on conservation laws, which 
must be gauge invariant. In order to reconcile this ap¬ 
parent contradiction let us explore more closely which 
quantities are gauge variant and which are the ones that 
must be conserved. 

With reference to Eq. <0. the measurable physical 
quantity that must be conserved by gauge transforma¬ 
tions is 771, which, in the gauge chosen, corresponds to 
the component of the angular momentum. However, 
it is crucial to note that is not (nor need be) gauge 
invariant. 

In fact the operator is defined as 


N 




X Pi 


(16) 


where p is the canonical linear momentum p = —iV. 
Although r is gauge invariant, p is gauge variant and 
therefore so is L^. In the following we wish to extend 


Eq. (101 so that the metric associated with the conserva¬ 


tion of 777 is indeed gauge invariant. 

We consider a system for which there exists at least 
one gauge such that = 0, with H the system 
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Hamiltonian. We name this the reference gauge and refer 
to its vector potential as A^e/ (r) and to its paramagnetic 
current density as (r). In this reference gauge the 

set {to} corresponds to the eigenvalues of and both 
equalities in the relation (11) hold. The set {to} is then 
a constant of motion and in this gauge it represents the 
z component of the angular momentum. 

We now focus on the generic gauge corresponding to a 
generic vector potential A (r). In this generic gauge, the 
first equality of Eq. (Ill holds, but the second equality 


holds only if is a constant of motion in this gauge. 
Here we consider the quantity 


jp W = jp W -pWVxre/ 

and the operator 

N 


L^ = ^[rx (p-VXre/)]^, 


(17) 


(18) 


where Vxref is defined by A = A^ef — We note 

that jp (r) is gauge invariant, as, from Eq. (|l4|. 


jp (r) = jp,.e/ (r) 

always. It follows that 


X jp (r) 


dr = TO 


(19) 


( 20 ) 


independently of the gauge. Furthermore, 

by using the definition 0 . where p (r) = 

/ |'0 (rr 2 ... rjv)|^ dr 2 ... drjv, and the first equality 
of Eq. (11 1 , which holds regardless of whether or not 


is a constant of motion, we obtain 


J r X jp (r) ^ dr = (V' I Tz I - y [r X p (r) Vxref]^ dr, 


= (V'lTz IV') - (^1 (r X VXre/)^ IV'), 
= ( 21 ) 


This demonstrates that Eq. (18 1 defines the operator as¬ 
sociated to the conservation law (201 independently of the 
gauge. In particular, comparison of Eqs. (20) and (21) 


shows that indeed is the operator whose eigenvalues 
are {to} independently of the gauge. [ 20 ] . 

Here reduces to in the reference gauge and in 
all gauges where is a constant of motion, as should be 
expected. This is because the limited set of gauges for 
which [Lz,H] — 0 holds, is the same within which both 
Lz and [r x jp (r)]^ are unaffected by gauge transforma¬ 
tions. These gauges correspond to vector potentials of 
the form 


A{r) = [xa + yl3,ya- ( 22 ) 

where a, j3, and 7 are all arbitrary functions of 
{x^ +y'^,z). These vector potentials are linked by gauge 


transformations of the form x ( 2 ;^ -|- , z). Demonstra¬ 

tion of these statements is quite long and the details are 
given in Appendix]^ Using the conservation law (20), 
we derive the metric, 



jpi (r) - jp 



2 


for the gauge invariant current density, jp. 


dr, 

(23) 


IV. METRIC ANALYSIS OF GROUND STATES 
FOR VARYING MAGNETIC FIELDS 

In the presence of a magnetic field the metric spaces for 
ground-state wavefunctions, particle densities and para¬ 
magnetic current densities are characterised by a “band 
structure” [ 8 ]. This is significant as identification and 
characterisation of ground-state properties is very impor¬ 
tant in several contexts but far from obvious. This “band 
structure” originates from the conservation law involving 
the paramagnetic current density jp. In jp metric space, 
when considering variations in the scalar potential [5] , the 
“band structure” is formed by spherical segments of al¬ 
lowed and forbidden distances on the concentric spheres, 
at least for the systems analysed. The specific arc length 
of these segments varies depending on the radius \m\ of 
the sphere. Here we wish to investigate how this “band 
structure” responds to changes in the magnetic held. 


A. Model systems 


We focus on two atomiclike model systems with uni¬ 
form, time-independent magnetic helds B = WcCZ ap¬ 
plied, where c is the speed of light. Both systems consist 
of two electrons in harmonic conhnement but with dif¬ 
ferent electron-electron interactions. One system, known 
as the magnetic Hooke’s Atom (HA), has a two-body 
Coulomb interaction lilll], whereas in the other elec¬ 
trons interact via an inverse square interaction (ISI), the 
relative strength of which can be varied through an in¬ 
teraction parameter a [22]. The Hamiltonians for these 
systems are 


Hha = ^ 


;A(u) 


L 2 2 
-0Jr,r, 




and 




n 2 


p* -b -A(r,) 
c 




1 

|r2 - ril’ 
(24) 


2 • 


(ri - r2) 
(25) 

Here, A = ^ (B X r) in the symmetric gauge, which is 
of the form of Eq. (22). Following Ref. [T8| , we neglect 


spin terms in the Hamiltonians to concentrate on the 
features of the orbital currents. For the ISI system, we 
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FIG. 2. (Color online) Energy plotted against the cyclotron 
frequency for several values of m for (a) Hooke’s Atom and (b) 
the ISI system. The confinement frequency and interaction 
strength are held constant. Arrows indicate where the value 
of m for the ground-state changes. 


can solve the time-independent Schrodinger equation ex¬ 
actly for all frequencies and values of m and a. However, 
for Hooke’s Atom, analytical solutions only exist for a 
discrete set of frequencies. In order to give us freedom 
over the frequencies we choose, we solve the Schrodinger 
equation with the method in Ref. [23) . which allows us 
to numerically determine the solution with high precision 
for all frequencies. 

We generate families of ground states by varying the 
magnetic field via the cyclotron frequency Wc, while hold¬ 
ing the confinement frequency, Wq, and all other param¬ 
eters in the Hamiltonian constant. For each Wc value we 
calculate the wavefunction, particle density, and param¬ 
agnetic current density. Within each family, one value 
of Wc (and hence m) is selected as a reference and 

niref respectively), with the appropriate metrics used to 
find the distance between the physical functions at the 
reference and all of the other states in the family. We 
choose the reference so that most of the available distance 
range is explored for both increasing and decreasing uJc- 

Figurej^shows that, for both of our systems, the value 


of m for which the energy is lowest decreases from zero 
through the negative integers as ojc increases. Conse¬ 
quently, when studying ground states, we must consider 
states on different spheres in the paramagnetic current 
density metric space. We also note that there are “tran¬ 
sition frequencies”, i.e., values of Wc where the energy 
is equal for two consecutive values of m. These are the 
crossings of the energy curves in Fig. Therefore, when 
varying uJc it is necessary to change the value of m at the 
“transition frequencies” in order to continue analysing 
ground states. Additionally, when m = 0, jp(r) = 0 
for all r. Hence, we take only negative values of m to 
ensure we consider ground states with nonzero paramag¬ 
netic current densities. 


B. Ground states’ band structure and relevance to 
current density functional theory 

An important research area where properties of ground 
states are central is DFT. This theory has produced 
widely used tools for realistic calculations of proper¬ 
ties of many-body systems, such as band structures of 
metals and semiconductors, crystal structures of solids, 
and characterisation of nanostructures [211 US] . DFT is 
founded upon the Hohenberg-Kohn (HK) theorem pS] . 
which states that there is a one-to-one mapping between 
ground-state wavefunctions and ground state particle 
densities. There are various forms of DFT that extend 
the application of the theory to a greater range of sys¬ 
tems. CDFT is the extension of standard DFT to include 
systems subject to external magnetic fields [TH|. There 
is a HK-like theorem at the core of CDFT that states 
that a one-to-one mapping exists between the ground- 
state wavefunction tp and, taken together, the particle 
density p (r) and the paramagnetic current density jp (r) 
(CDFT-HK theorem) [TB]. This additional complexity of 
the CDFT-HK theorem with respect to the original HK 
theorem is due to the fact that systems with magnetic 
fields are characterised not only by a scalar potential 
(the external potential), but also by the vector poten¬ 
tial connected to the magnetic field [18]. In Ref. [8] we 
started examining the CDFT-HK mapping by looking at 
the effect of varying the scalar potential, i.e., the external 
confining potential; here we wish to complete this anal¬ 
ysis by looking at the effect on the mapping of varying 
the vector potential, i.e., the magnetic field. 

We start by comparing the distances between wave- 
functions, their related particle densities, and their re¬ 
lated paramagnetic current densities. Figure [^ shows 
plots of the relationships between the various distances 
considered, with each point referring to a particular value 
of oJc- Let us consider first the plots of particle density 
distance against wavefunction distance [Figs.j^a) -[^d)]. 
As observed in Ref. [8] , metric space regions correspond¬ 
ing to ground states present a “band structure”, where 
points associated with the same value of |m| are grouped 
into distinct segments, i.e., bands. However, in contrast 
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FIG. 3. Plots of distances for Hooke’s Atom with reference state cjq = 0.5, cac = 5.238, and m^e/ = —5 (top two rows) and for 
the ISI system with reference state ujq = 0.6,tJc = 5.36,Of = 5, and mrej = —10 (bottom two rows): (a) - (d) particle density 
distance against wavefunction distance, (e) - (h) paramagnetic current density distance against wavefunction distance, and (i) 
- (1) paramagnetic current density distance against particle density distance. The reference frequency is taken halfway between 
the two “transition frequencies” related to niref- 


to the band structure observed in Ref. [5] [sketched in 
Fig. Qa)], when varying the vector potential we obtain a 
series of “overlapping bands”, where the minimum wave- 
function and minimum particle density distances for one 
value of |to| are smaller than the maximum distances for 
the previous value of [mj. This implies that there is an 
overlap between the projections of the bands on the met¬ 
ric space sphere representing the densities, as sketched 
in Fig. I^b) (similarly for the projection on the sphere 
representing the wavefunctions). Though overlapping, 
this band structure still results in discontinuities in the 
relationship between Dp and when the value of m 


changes. Unlike when varying loq [H], by varying the mag¬ 
netic field we do not observe any forbidden distances, so 
we cannot identify forbidden regions for ground states by 
considering the particle density and wavefunction metric 
spaces alone. In the range of distances explored here, 
nearby wavefunctions are mapped onto nearby particle 
densities and distant wavefunctions are mapped onto dis¬ 
tant particle densities. However, in contrast to Ref. [5], 
the mapping is only piecewise linear: When acting on 
the vector potential, as Wc is swept through each tran¬ 
sition frequency, ground states and their particle densi¬ 
ties abruptly revert to be closer to the reference state. 
























7 


(a) (b) 



N=2 N=2 

FIG. 4. (Color online) Sketches of “band structures” consist¬ 
ing of (a) “bands” and “gaps” and (b) “overlapping bands” 
in particle density metric space for three consecutive bands, 
where a different patterning corresponds to a different value 
of m. The reference state is at the north pole. 


while an almost linear mapping is maintained within two 
consecutive transition frequencies. The segments created 
in this way do not overlap, as, at each transition fre¬ 
quency, the ball related to the particle density and cen¬ 
tered at the reference density shrinks proportionally more 
than the corresponding ball related to the wavefunction. 
Also, in contrast with Ref. [8] , the two families of ground 
states corresponding to \m\ < luire/l and |m| > \mref\ 
describe distinct paths in metric space [e.g., compare 
Figs. [3|;a) and|3];b)], with the size of the bands greater 
for |m[ < \mref\ compared to \m\ > \mref\- For all 
of these reasons the CDFT-HK mapping between wave- 
functions and related particle densities acquires added 
complexity when varying the vector potential compared 
to varying the scalar potential [compare Figs. 2(a) and 
2(b) in Ref. [5] with Figs, [^a) -j^d)]. 

In Figs. I^e) - [^h) we consider paramagnetic current 
density distance against wavefunction distance. Here we 
find once more an overlapping band structure for wave- 
function distances; however a band structure with re¬ 
gions of allowed (bands) and forbidden (gaps) distances 
is observed for paramagnetic current density distances. 
In contrast with the one sketched in Fig. |^a), in this 
case each band resides on a different sphere according to 
the value of |to| (the radius of the sphere). Transition 
frequencies are points of discontinuity for both paramag¬ 
netic current density and wavefunction distances. As for 
Figs. ia)-id), the curves for increasing and decreasing 
ujc (and hence \m\) do not overlap, with larger bands for 
small values of \m\. Finally Figs. ii) -il) present the 
plots of paramagnetic current density distance against 
particle density distance. These exhibit behavior similar 
to that in Figs. §e)-§h). 

The overlapping band structures observed in Fig. 
demonstrate that mappings between some of the dis¬ 
tances considered here are multivalued. This multival¬ 
ued mapping does not represent a contradiction of the 
CDFT-HK theorem as it is entirely possible to have dis¬ 
tinct functions at the same distance away from a refer¬ 


ence. In particular, in terms of the “onion-shell” geom¬ 
etry, all states situated at the same polar angle and on 
the same sphere will have the same distance from the 
reference state. 

In Fig. the wavefunction and paramagnetic current 
density distances are plotted against Wc for both systems, 
enabling the band structures for individual functions to 
be analysed. We note that, as observed in Fig. there 
is a decrease in the wavefunction distance at transitions 
[Figs.a) and[^b)], but an increase in the paramagnetic 
current density distance [Figs. ie) and[^d)]. These fea¬ 
tures give rise to overlapping-band and band-gap struc¬ 
tures, respectively. The other major feature is that, when 
varying Wc, there is nonmonotonic behavior within bands 
corresponding to values of m close to m^e/ (see insets). 
For both wavefunctions and paramagnetic current den¬ 
sities, we observe that immediately after each transition 
frequency, the distances initially decrease to a minimum 
for that particular band before increasing to the max¬ 
imum for the band. This occurs at the transition fre¬ 
quency to the next band. This behavior is more pro¬ 
nounced for wavefunctions than for paramagnetic current 
densities. As stated, the nonmonotonicity is not in con¬ 
tradiction with the HK-like theorem of CDFT, but shows 
a richer behavior with respect to what was observed in 
Ref. [5] when varying the scalar (confining) potential. 

We point out that the band structure in metric space 
for paramagnetic current density is fundamentally differ¬ 
ent from the ones for particle density and wavefunction, 
as the former develops on different spheres, one band for 
each sphere, while the latter are within a single sphere 
where they may display overlapping-band or band-gap 
structures (see Fig. |^. All these band structures origi¬ 
nate from the conservation law characterising the para¬ 
magnetic current density and the features of the metric 
spaces for wavefunctions and particle densities are a di¬ 
rect consequence of the mapping of jp (r) onto ip (r) and 
onto p(r). In this sense the band structure features of 
the metric spaces for wavefunctions and particle densities 
could be seen merely as the projections done by these 
mappings of the band structure characterising the para¬ 
magnetic current density. 

Finally, we wish to concentrate on the implications of 
our findings for CDFT. CDFT requires that both p and jp 
are taken together to ensure a one-to-one mapping to the 
wavefunction. The metric analysis allows us to provide 
evidence for an important aspect of this mapping, that 
is, to understand when the inclusion of jp in the mapping 
becomes really crucial for the one-to-one correspondence 
to hold. 

We present in Fig. I^the ratio /Dp against for 

both Hooke’s atom and the ISI system. From the data 
it is immediately clear that, in metric space, to a good 
level of approximation, = const x Dp as long as 

m = niref- This constant is the same for lOc > 
and ojc < These findings suggest that, at least 

for the systems at hand, as long as we remain on the 
same sphere in the paramagnetic current density metric 


















FIG. 5. For Hooke’s Atom (top) and the ISI system (bottom), (a) and (b) wavefunction distance and (c) and (d) paramagnetic 
current density distance are plotted against tJc- The behavior around the reference frequency is shown in each inset. The 
reference states are ujo = 0.5, = 5.238, and m^e/ = —5 for Hooke’s Atom and cjo = 0.6,cjc,.^f = 5.36, a = 5, and 

rriref = —10 for the ISI system. 


space, jp and p carry very similar information and the role 
of jp in the core mapping of CDFT is secondary. The 
situation becomes very different for ground states with 
m yf TTiref ■ In this case the ratio /Dp is far from 
constant and Fig. clearly shows that the information 
contents of jp and p are both necessary to define the state. 
Similar results are obtained when keeping the magnetic 
field fixed but varying the confinement ujq of the systems 
(not shown). 

The characterisation of this difference in the role of jp 
and p in the CDFT core mapping constitutes one of the 
main results of the paper. To support it, we will analyse 
in the next section the behavior of states where m is kept 
equal to niref at all values of uic- 

V. EXCITED STATES 

Although an understanding of the ground state is im¬ 
portant for studying systems subject to magnetic fields, it 
is often necessary to go beyond ground states, for exam¬ 
ple, when studying rapidly varying fields or spintronic de¬ 


vices that operate with excited states. With the metrics 
at hand, we investigate excited states and consider dis¬ 
tances between families of states corresponding to fixed 
values of m m- For each value of m we will construct a 
family of states by varying uic (wo and a kept constant) 
and calculating the corresponding wavefunctions, parti¬ 
cle densities, and paramagnetic current densities. As 
for ground states, we choose m < 0. With respect to 
Fig.i this corresponds to following single energy curves 
smoothly, i.e., without switching to a different curve at 
crossings, as done instead for the ground state case. Each 
family of states will then lie on a particular sphere in 
the paramagnetic current density metric space. As the 
states considered are not necessarily ground states (see 
Fig.[|) , there is no one-to-one mapping between the wave- 
function and particle and paramagnetic current densities, 
but, these being fundamental quantities that characterise 
the system, we will still explore their relationships. Ad¬ 
ditionally, the study of these quantities allows us to cor¬ 
roborate the findings related to Fig. 

Figure shows the relationship between each pair of 
distances for six different values of m. For any pair of 
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FIG. 6. Plots of the ratio of paramagnetic current den¬ 
sity distance to particle density distance against wavefunc- 
tion distance for (a) Hooke’s Atom with reference state cuo = 
0.5,cJc^^f = 5.238, and m^e/ = —5, and (b) the ISI sys¬ 
tem with reference state ojo = = 5.36, q = 5, and 

mref = — 10 . 


distances here discussed, we find a monotonic relation¬ 
ship that is linear in the short- to intermediate-distance 
regime, before one of the two functions rises more sharply 
to its maximum (see also Fig.[^. The mapping between 
the physical functions is such that nearby functions a 
(e.g., the wavefunctions) are mapped onto nearby func¬ 
tions b (e.g., the paramagnetic current densities) and dis¬ 
tant functions a are mapped onto distant functions b. 
Crucially, as opposed to ground states, distances do not 
form any kind of metric space “band structure”, confirm¬ 
ing the origin of band structures as the changes in m. 

Looking at wavefunction distances against particle 
density distances in Figs. [^a) and [^b), and contrast¬ 
ing with Figs. [^a)-[^d), we observe that the curves for 
increasing and decreasing uic and all values of |m| collapse 
onto one another. This hints at a universal behavior of 
the mapping between particle density and wavefunction 
when all the physical quantities describing the system 
remain on the same sphere in the related metric space 


while a physical parameter is smoothly changed. 

When considering paramagnetic current density dis¬ 
tance against wavefunction distance in Figs. [^c) and 
Hd), although the curves for increasing and decreasing 
LOc collapse onto one another, the curves for different val¬ 
ues of m are distinct, particularly when \m\ is low. For 
lower values of \m\ the linear region extends across a 
larger range of distances. There is also a relatively small 
increase in the gradient at greater distances for low |m|. 
The curves in Figs.j^c) and[^d) all start and end at the 
same points. With the rescaling for used in Fig. [Tl 
the curves tend to a limiting curve with increasing vff 
ues of \m\. In Fig. we show the relationship between 
wavefunction distance and paramagnetic current density 
distance for the ISI system without rescaling . Here, 
the curves for each value of |m| intersect only at the ori¬ 
gin, and each has a unique maximum of 2 \m\ for the 
paramagnetic current density distance. We observe that 
the gradient of the initial linear region increases with \m\. 
Figure shows, for the ISI system, that the gradient in 
this region increases linearly with |m|, « fc |to| 

with 0 ^ fc ^ I, and is approximately equal for both de¬ 
creasing and increasing ujc- Similar results are obtained 
for Hooke’s Atom (not shown). These results imply that 
when rescaled as in Fig. the initial slope of the curves 
will always be below 45°, a result also observed in Ref. 
for the case in which different spheres in the wavefunction 
metric space geometry were considered. 

When considering paramagnetic current density dis¬ 
tance against particle density distance [Figs. [^e) and 
we see that, as for vs with the rescaling 

m Fig. [^there are distinct curves for each value of m that 
converge onto a single curve as |to| increases. As opposed 
to vs D^p, the extent of the linear behavior of these 
curves is increasing as jmj increases. 

The behavior of the curves observed in Fig. [^reflects 
the “onion-shell” geometry. For wavefunctions and par¬ 
ticle densities the sphere radius is associated with the 
number of particles in the system, which is fixed for the 
systems considered. Thus, regardless of the value of jmj, 
wavefunctions and particle densities always lie on the 
same sphere in their respective metric spaces. The fact 
that the related curves still superimpose for changing |to| 
seems to imply that the value of jmj has no relevant ef¬ 
fect on the curves representing the relative change of '0 
and p for changing parameters, at least as long as they 
remain on the same sphere. In paramagnetic current den¬ 
sity metric space, the spheres’ radii are related to jmj, so 
paramagnetic current densities are on the surface of dif¬ 
ferent spheres each time we consider a different value of 
jmj. As a result we see that the curves’ shape is affected 
and they do not collapse onto each other. A similar uni¬ 
versal behavior within each sphere and, by contrast, the 
breaking of this universality when different spheres were 
considered, was also observed in Ref. [3, where different 
values of N, and hence different spheres, for both wave- 
functions and particle densities were considered. This 
seems to suggest that different behavior for the mappings 
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FIG. 7. Plots of: (a) and (b) particle density distance against wavefunction distance, (c) and (d) paramagnetic current density 
distance against wavefunction distance, and (e) and (f) paramagnetic current density distance against particle density distance 
for m = —1, —2, —3, —8, —9, —10. The reference states, for each value of m, are: loq = 0.1 and = 30.0 for Hooke’s Atom 
(top) and ujo ~ 0.1, = 5.0, and a = 5 for the ISI system (bottom). Closed symbols represent decreasing lUc and open 

symbols represent increasing ujc- 


should be expected when curves on different spheres in 
the metric spaces are involved. 

Finally, Fig. combines all distances for each sys¬ 
tem in a single plot. Importantly, this figure shows that 
for small to medium wavefunction distances /Dp ^ 

constant, where the constant depends on |m|, so this 
ratio is, to a good approximation, independent over vari- 



FIG. 8. Plot of paramagnetic current density distance against 
wavefunction distance for m = —1,—2,—3,—8,—9,—10 for 
the ISI system. We take the state with too = 0.1,u)c,.^f = 5.0, 
and a = 5 as a reference for each value of m and consider 
distances across the surface of each individual sphere. 


ations of the wavefunction for relatively close wavefunc- 
tions. In this respect, for relatively close wavefunctions 
this suggests that the mappings between current density 
and wavefunction and between particle density and wave- 
function are very similar, as long as the family of states 
follows the evolution of the same energy eigenstate as 



FIG. 9. Plot of the ratio of paramagnetic current density 
distance to wavefunction distance against |m| for the ISI sys¬ 
tem. The reference for each value of m is cuq = 0.1, = 

5.0, 0 = 5, and the gradient is taken at ljc = 4.5 for decreas¬ 
ing frequencies and Uc = 6.0 for increasing frequencies, i.e., 
the frequencies corresponding to the closest points to 
for both decreasing and increasing cUc in Fig. 
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FIG. 10. Plots of the ratio of paramagnetic current den¬ 
sity distance to particle density distance against wavefunc- 
tion distance for (a) Hooke’s Atom, with reference state 
cjo ~ = 30.0, and (b) the ISI system, with refer¬ 

ence state uiQ = 0.1,a;c^,y: = 5.0, and a = 5. Closed symbols 
represent decreasing ujc and open symbols represent increas¬ 
ing LOc- 

driven by the varying parameter (see Fig.j^. 

VI. CONCLUSION 

The metric space approach to quantum mechanics has 
enabled us to illustrate the role of the vector potential in 
systems subject to external magnetic fields, with partic¬ 
ular reference to the fundamental concepts of CDFT. Im¬ 
portantly, we have also furthered the theoretical frame¬ 
work of the metric space approach to quantum mechan¬ 
ics by discussing the key point of gauge invariance for 
the natural metrics proposed and in particular demon¬ 
strating the gauge invariant metric for the paramagnetic 
current density, which is not a gauge invariant quantity 
per se. 

The presence of the vector potential in the Hamiltonian 
leads to the inclusion of the paramagnetic current density 


in the core mapping of CDFT. By considering the met¬ 
ric for the paramagnetic current density together with 
that for the particle density, we were able to investigate 
the relative contribution to this core mapping from each 
of these two quantities, for the systems at hand. When 
m is held constant, and paramagnetic current and parti¬ 
cle densities belong to the same metric space sphere as 
their reference state, we observed that the ratio /Dp 
is approximately constant, suggesting that p and jp con¬ 
tribute similar information. However, this simple relation 
dramatically breaks down when considering states with 
m yf rriref and hence states spanning different spheres 
in paramagnetic current density metric space. This sug¬ 
gests that the presence of jp in the core mapping of CDFT 
becomes crucial in this case. 

By varying the vector potential, we uncovered different 
aspects of the “band structure” in ground-state metric 
spaces, in particular the presence of overlapping bands, 
which enriches the band-gap structures already observed 
when varying the scalar potential. Our analysis suggests 
that, in general, the presence of bands in metric space can 
be expected when considering a family of states for which 
one of the fundamental physical functions spans more 
than one sphere in its metric space. For ground states, 
the onset of the band structure is the signature of energy 
levels’ crossings obtained by varying a parameter in the 
Hamiltonian (the magnetic field in the present case). 

We also applied the metric space approach to quan¬ 
tum mechanics beyond ground states. When consider¬ 
ing families of states characterised by fixed values of to, 
it was found that the mappings between wavefunctions, 
particle densities, and paramagnetic current densities are 
monotonic and almost linear, without band structures, 
confirming that each band is characterised by a specific 
value of TO. The curves Dp versus Dp superimpose for 
all values of to, but not so when D-, is involved. This 
is consistent with the fact that a different to represents 
different spheres in the “onion-shell” geometry related to 
jp- 

Finally, when considering the ratio Dj^^ /Dp for these 
“fixed-TO” families, the relationship /Dp « const 

was observed to persist up to intermediate distances, and 
for all of the values of to that were explored. At least for 
the systems at hand, this suggests that, within the same 
sphere and up to quite different states, particle density 
and wavefunction still suffice to contribute most of the 
information on the physical system. 
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Appendix A: Gauge invariance for the wavefunction 
metric 

Gauge transformations affect wavefunctions by intro¬ 
ducing a constant global phase factor [see Eq. (HU)]- 
Wavefunctions differing by this phase factor describe the 
same physics; in fact, the solutions of the Schrddinger 
equation are only defined up to a global phase factor. 
To have physically meaningful metrics, it is therefore im¬ 
portant to define equivalence classes such that the metric 
assigns zero distance to wavefunctions differing only by 
a global phase factor. 

An equivalence class for an element a; S A is defined 
as pH] 


With this semimetric space ({'0},!?^), we can recover 
a metric space in a natural way, by “gluing” equivalent 
elements to form a set of equivalence classes. By con¬ 
sidering the set of equivalence classes, rather than the 
set of all wavefunctions, all wavefunctions differing only 
by a global phase factor are identified with one another. 
Thus, for wavefunctions, the set of equivalent wavefunc¬ 
tions with is a metric space, with the metric defined 
between each of the equivalence c lasse s, as required [SH] . 
The metric D^p defined from Eq. (A 61 is the same as the 
metric defined in Refs. liiHni arid can be expressed in 
the form of Eq. ® . 


Appendix B: Determining the gauges where is a 
constant of motion 


[x] = {x' G A : X ~ x'} , (Al) 

where ~ is the equivalence relation. Each element of the 
set A belongs to a single equivalence class [28] . 

In order to account for an equivalence relation between 
elements, x ^ x', we follow a general procedure for intro¬ 
ducing equivalence relations into a metric space {X,D). 
We define the function [29] 

DR{x,y) = inf = x, = ?/,/c G n| , 

(A2) 

where the infimum is taken over all choices of {pi} , and 
{qi} such that Pi+i ~ Qi- This implies that if x ^ y. 
Dr (x, y) = D (x, x) + D {y, y) = 0 even if D (x, y) ^ 0 
[29] . This function is a semimetric (or pseudometric) on 
the set A, known as the quotient semimetric. A semi¬ 
metric is a distance function that obeys all of the axioms 
of a metric except that it allows zero distance between 
nonidentical elements as well as identical ones. 

For wavefunctions, the metric derived from the conser¬ 
vation law before accounting for the equivalence of wave- 
functions differing by a global phase is 0(30] 


In order to be a constant of motion, the 2 component of 
the angular momentum = —i [r x V]^ must commute 
with the Hamiltonian. Given that a vector potential is 
present, we consider the Pauli Hamiltonian 


H = —V-" 
2 


ihe 

2me 


A-V-f V-A- 


2m, 


with V (r) such that V (r), 


= 0 . 


-H(r), 
(HI) 

The Hamilto¬ 


nian (HI) does not necessarily commute with for a 
particular A(r), because is gauge varia nt. For in¬ 
stance, Lz commutes with the Hamiltonian (Bl| in the 
symmetric gauge A = [y, —x, 0] and does not commute 
with it in the Landau gauge A = [0, —x,0]. We wish 
to determine the general set of vector potentials where 


H,Lz 


= 0. ET] 


1. Simplifying the commutator 

The commutator we wish to evaluate is 


Dp, 




(A3) 


For this in general we have that Dp, (i/), ^ 0. If in 

Eq. ( [A2| ) we take fc = 2, we find 

Dp, = inf ^Dp, + Dp, (V’ 2 ,'*/' 2 )| (A4) 

= inf (■0i,'0')} ■ (AS) 

where ip' = D'^' 4’2 ~ V '2 and we have used the positivity 
axiom of the metric. The choice of ip' that will minimise 
the value of the semimetric is determined by the phase 
factor, hence 

Dp,{ipi,ip 2 ) = miii^Dp, {ipi,e"^ip 2 )'^ (A6) 



+ 4(a.v« + v.(ax|)-v.(a,|^ 

r) f) Filh 

- x^X ■ {AiP) + • {AiP) + A^x"^ 

- {A^ip) + y^ (A'^ip ), 

OX oy ^ ^ ox ^ ^ 

(B2) 


where we have used that 


^V 2 + H(r),L, 


= 0. We 


H,L, 


= 0, and then solve 


wish to impose the condition 

the commutator to obtain the vector potential A(r). 
After performing the vector operations and simplifying. 




















13 


Eq. (B2) reduces to 


oy ax j \ax ay 


dip dAy dip dA^ 


dy dy dz dy 


2y 


— xip 


dip dAx dip dAy dip dA^ 
d (dA^ 


dy \ dx 
d (dA^, 


dy dx dz dx 
dAy d^ 
dy dz 


dA,, dA, 


dx V dx 


dy dz 


, ^242 , dA^ 


(B3) 


In order to progress with the solution of this equation, 
we first consider the case where ip, and ^ are 

all independent of each other. This choice allows us to 


decompose Eq. (B3) into a set of simultaneous equations, 


which we can then solve. The solution of these equations 
will provide properties of the general set of vector poten¬ 


tials where 




= 0. Using these properties, we will 


then solve Eq. (B31 for A (r) using a general wavefunc- 
tion. 

With our choice of trial wavefunction, we write the set 
of simultaneous equations 


— X 


d fdAx dAy 


dy \ dx dy 


dA, 

dz 


— X- 


+ 2/7^ 


d (dA^ 


dAy dA^ 


dx V dx 
dA 


dy 
dAx 

Ax - x^A^ + y^r^ = 0 , 


dz 


Ay -fa; 

dy 
dA„ 


dy 


dx 


dA,, dA,, 

y— - X —— = U. 


dx 


dy 


dA^ 

dy 

-.^=0, (B4) 

(B5) 
(B6) 
(B7) 


We concentrate first on Eqs. (B5)-(B7|, a set of three 
equations for the three unknowns Ax, Ay, and A^. 
Firstly, we consider Eq. ( |B7| ): In order to solve this par¬ 
tial differential equation (PDE), we use the method of 
characteristics [32]. 

The method of c haracteristics requires the visual¬ 
isation of Eq. (B7) in four-dimensional coordinates 
{x,y,z,u). By considering the solution surface u = 
A, {x, y, z), we can write 

A, {x,y,z) - u = 0. 

For any surface, S, a normal vector to the surface is given 


by VS. Thus, the vector 


dA, dA, dA, 
dx ’’ dy ’’ dz 


-1 


is normal 


to the solution surface. We now write the PDE (B7) as 
a scalar product 


[y, -X, 0, 0] 


dA, dA, dA, 
dx ’ dy ^ dz 




= 0 . 


Since the scalar product of these two vectors is zero, 
they must be orthogonal. Given also that the vector 


-1 


is normal to the surface, this tells 


dA, dA, dA, 
dx dy ’’ dz ^ 

us that the vector field [y, —x, 0,0] is tangent to the sur¬ 
face at every point, providing a geometrical interpreta¬ 
tion of the PDE. Thus, any curve within the surface 
Az{x,y,z) — u = 0 that has the vector [y, —a;,0,0] as 
a tangent at every point must lie entirely within the sur¬ 
face. Such curves are called characteristic curves [32] . 
Any curve can be described by a parameter t and the tan¬ 
gent of such a curve r (t) is given by the derivative with 
respect to this parameter r' (t). Therefore, the tangent 
of a characteristic curve r (t) = [a; (t) ,y{t),z (t), A, (t)] 
is given by the vector 


r' (t) = 


dx 
dt ’ 


dy 

dt 


dz dA, 
dt ’ dt 


This vector is therefore proportional to the tangent vec¬ 
tor [y, —a;,0,0] for this characteristic curve, allowing us 
to construct the equations 

dx 
dt 
dy 
dt 
dz 
dt 
dA, 
dt 


= y, 

(B8) 

= -X, 

(B9) 

= 0, 

(BIO) 

= 0. 

(Bll) 


These are the characteristic equations of the PDE (B7). 
Solving this set of ordinary differential equations 


(ODEs) yields the solution of the original PDE (B7|, 
since 

dA, dx dA, dy dA, dz dA, 

dt dt dx dt dy dt dz 

dA, dA. 


dx 


— X 


= 0 . 


By eliminating the parameter t in Eqs. (B8)-(B11), we 


can reduce the set of ODEs to three equations 


dy X 

dx y ’ 

(B12) 

dz 

(B13) 


dA, 

(B14) 

-r^ = 0. 

dx 


We now note that the constant of integration in 


Eq. (B14) has a functional dependence on the solutions 
to Eqs. (B12) and (B13). This is because the ODEs are 


solved along characteristic curves: The constants of in¬ 
tegration are constant along a particular characteristic, 
but can vary between characteristics. The solutions to 


Eqs. (B12) and (B13) are 


2 I 2 
X + y = a 


z = b. 


(B15) 
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respectively, where a and b are the constants of integra¬ 
tion. Thus, the solution for A, is 


Az = j {x'^ + , z ), 

where 7 is an arbitrary function. 

2. Solving the simultaneous equations 


(B16) 


of the characteristic equation, we have complete freedom 
in the choice for rj, provided that it is not the same as 
Given that we know that the symmetric gauge satis¬ 


fies the commutator 


H,Lz 


= 0 and, specihcally, that 


Ax = y is a solution for Eq. (B19l, we choose 77 = y. By 


using the chain rule, we find the derivatives in Eq. (B19l, 


We will now solve Eqs. (B5) and (B 6 ) simultaneously. 


First, we differentiate Eq. (B5) wit^ 
and y, which gives 


dAx „ dA 
= - 

ox 
dAx 


1 respect to both x 


dy 


= 2y 


dAx 


dAx 


dr] 


dx 

dAy 

dy 


-I- X 


dAx 

dy 

d^Ax 


dy'^ 


d^Ax 
dxdy 
_ dAx 
dx 




y 


d^Ax 
' dx^ 
d^Ax 
dxdy 


= 0 , 


respectively. We substitute these expressions for and 



d^Ax 


r^dAx 
+ 2 , 

(B17) 

dx"^ 

de 

d^ ’ 

d'^Ax 

, 2d^Ax 

. 9^ Ax 

(B18) 

dy"^ 

= 4u^ 
y g^2 

^'^^d^dy 

d^Ax 

. d'^Ax 
= 4xv - r- 



dxdy 

y Q^2 

d^dy 


d'^Ax dA 

- - -I- 2 

dy^ 


d^ ’ 


dA^ 


into Eq. (|B 6 |) and obtain 


Ax ^ 

d^Ax dAx 

d^Ax 



^ dy"^ dx 

^^dxdy 

) 


f 

dAx d'^Ax 

d^Ax 

)=o. 


+ 2 / (- 

dy ^ dxdy 



29"^ Ax 

. d^Ax 

.d^Ax 

dAx 

dAx 

y dx^ 

oxoy oy^ 

^ dx 

^ dy 


Substituting into Eq. (B19l and simplifying, we get 

nd'^Ax dAx . 

x^^^-y^+Ax = Q 
drj^ dr] 

and completing the change of variables gives 


+ Ax 
(B19) 


We now have an equation containing only the unknown 
Ax that we can solve. 

We begin to solve this equation by using the method of 
characteristics. For second order PDFs, it is first neces¬ 
sary to determine the type of the PDE, either hyperbolic, 
parabolic, or elliptic. This is done by calculating the dis¬ 
criminant b^—4ac, where a, 6 , and c are the coefficients of 
^dxdy ’ respectively. This will then allow 

us to perform an appropriate change of variables from 
(x, y, z) to (^, y, z), where ^ and t] are the characteristics 
1831 . The discriminant is 


= 0 . 




d^Ax 

di]'^ 


V- 


dAx 

dr] 


■ Ax — 0 . 


(B22) 


The next step of the solution is to perform a reduction 
of order through the use of the known solution Ax = t] 
[321 ■ The reduction used is 


Ax = UT], 


dAx du 
dr] dr] 


d'^Ax _ d^u du 
dr]'^ ^ dr]'^ dy ’ 


which we substitute into Eq. (B22) to give 


6 ^ - 4ac = 4x^y^ - 4a;^y^ = 0. 


(B20) 


Therefore, the characteristic equation is parabolic and 
has one repeated solution, which we take for The 
characteristic equation is the ODE [33] 


9 , d'^u „ du ndu 

,, d'^u du 

We now make the substitution v = ^, 


solving for ^ gives 


— 2xy^ -I- = 0; 

dx 


(B21) 


This equation can now be solved separably, 

dy. 


(B23) 


[ -dv= [ 

\ W - 2 C 1 

1 V J 

[y(^- y 2 )J 


dx 


X 

y' 


Hence, from Eq. (B15) we know that the hrst character¬ 
istic is ^ = a = sP~+ y^. Since there is only one root 


We decompose the denominator through the use of par¬ 
tial fractions, giving 

fldv = - f-dy+ [ -^dy, 

J V Jr] J i-y^ 
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from which we get 


InII = —21n7y — 2 I? ~ ^^1 + In [a (^, z)], 
a (C, z) 


rj2 1^- ry2|2 


where a is an arbitrary function and we note that = 
> 0 , hence ^ is always positive. 


Now that we have a solution to Eq. (B23l, we must 


reverse our substitutions to get a solution for A^. First, 
we integrate v to get u, 


u = a{^,z) [ - - ^ - rdy- 

J r)^ yy _ „2'i 2 


jf — rf-y 
From standard integrals [31], we get 
a{i,z) 


(B24) 


u = -- 




+Hy,z)., 


where /3 is another arbitrary function. Next we write 


Ax = ur], obtaining 

Ax = a (^, z) -r^Y + vP (^, z) 


(B25) 


where we absorb the factor of — | into a. Finally, we 
substitute back from {^,r],z) to {x,y,z), 

Ax = a {x^ + z) {x^ + y^-y^Y+ yP (x^ + y^z), 
= xa{x‘^+y‘^,z)+yl3 {x'^+ y‘^,z) , (B26) 


to give the solution for Ax- We find Ay from Eq. (B5|, 


Ay =y 


dAx 


dAx 


^ ^ dx dy 

=y [a {x^ + z) + 2x^a' {x^ +y^,z) + 2xyfi' {x^ + y^, z)] 

- X \2xya' (a;^ + y^, z) + /3 (a;^ + y^, z) + 2y'^l3' (x^ + y^, z)] 
=ya {x^ + y^, z) - xP [x^ + y^, z) , 

giving us solutions for all three components. 

We will now verify that the solutions for AxTAy, and 
Az, satisfy the remaining equation (B4), 


d_ / dAx 9^ dA^ 

dy\dx dy dz 

d 


dA^ d fdAx dAy dAx\ dA^ 


= — X— 1 2 q! + 2 a:^a' + 2xyP' + 2y^a' — 2xyP' + ^ 
oy \ oz 


y-^ f 2 a + 2 a;^a' + 2xyP' + 2y^a' — 2xyP' + ^ 
dx V dz 


— X {Ax^yaa' + 2y/3^ + 4y^/3/3' + 2ya^ + Ay^aa' + Ax'^yPP' + 4y77') 

+ y {2x0^ + Ax^aa' + Axy'^PP' + Axy^aa' + 2a;/3^ + Ax^PP' + 4 x 77 ') 

= — Axya — Ax^ya" — Axya — Axy^a” — 2xy^- + Axya + Axya + 4a:^ya" + 4a;y^a" + 2xy^- = 0, 


where a prime denotes differentiation with respect to [x^ + y^). 


Therefore, the form of the vector potential required for 



0 is. 


A = [a:a (a;^ + y^, z) + yP (a;^ + y^, z) , ya {x^ + y^, z) - x/3 (x^ + y^, z) , 7 (x^ + y^ z)] . (B27) 


This form of the vector potential satisfies the condition 


H,L, 


= 0 for wavefunctions with 


d'lfj dip 
dx ’ dy 


and 


dip 

dz 


all 


independent of each other. Clearly, vector potentials that are not of the form of Eq. (B27) do not satisfy the condition. 
However, in order to ensure that vector potentials of this form satisfy this condition for an arbitrary wavefunction, we 
use the properties of the vector potentials to solve the original commutator Eq. (B3) for an arbitrary wavefunction. 
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This gives 


H,L, 


4> =2 {xa + y/3)^ -2 {ya - xjS)^ - 2x 
oy ox 


{2xya' + /3 + 2y'^l3') ^ + {'^ + - 2a:y^') ^ 


+ 2y 


(a + 2x^a' + 2xyP') ^ + (2xya' — /3 — 2x^/3') ^ + 2 x 7 '^ 
' ' ax ^ ' ay oz 


— xip-^ f 2a + 2a;^a' + 2xyP' + 2y^a — 2xyP' + ^ ) 

ay \ azj 


d 


dx 


dj 


+ ytlj— 2a + 2a;^a' + 2xyl3' + 2y^a' — 2xyP’ + — 


dz 


— xtjj {Ax^yaa' + 2y^^ + Ay^jSP' + 2ya^ + Ay^aa' + Ax‘^y/3f3' + Ayj^') 

+ j/V' (2xa^ + Ax^aa' + Axy'^P/3' + Axy^aa' + 2xj3‘^ + Ax^jSfi' + 4 x 77 ') . 


This simplifies to 


H,L, 


d 


dj 

dz 

d ( Q'~s 

+ yi>-^ ( 2a + 2x^a' + 2y‘^a' + — 


lb = — x'tpTT 2 a + 2 x a + 2y a + 
dy 


r 


gauge transformations, /3 must be constant in a gauge 
transformation. 

A gauge transformation is given by A' = A + Vx and 
takes the form 


dx 


dz 


= — Axya' — Ax^ya" — Axya' — Axy^a" — 2xy^^ 
+ Axya' + Axya' + Ax^ya" + Axy^a" + 2xy 


dj^ 

dz 

di 


Vx = A' - A, 

= [xAa + j/A/3, j/Aa — xA/3, A 7 ] 

= [xAa, j/Aa, A 7 ], (B28) 


= 0 . 


r)y 

using that A/3 must be zero. We obtain x by integrating 
each of the components of the vector (|B28[) : 


Thus, the vector potentials of the form (B27l fulfill the 
condition =0 and in these gauges is a con¬ 

stant of motion. 


3. Gauge transformations between vector 
potentials for which Lz is a constant of motion 


xAadx = 


X(x^ +y^,z) 


j A-ydz = V {x^ -|- j/^, z) . 


,,, -1 . r ; • 1 , . Clearly then, the scalar field x must be a function of the 

We now consider a ga uge transformation between two . ( 2 2 \ 

ancpaa r^f fnrm ^Tt97^ A T^nfonl-icil r»f fViia fnrm ^ X V ^ 

Finally, we demonstrate that [r x jp]^ is unchan ged b y 
gauge transformations between gauges of the form (B27|. 


gauges of the form (B27). A vector potential of this form 
gives the magnetic field 


B = V X A 

= V X [xa -I- yd, ya — xj3, 7 ] 

„ / da dp da d/3 ^ , 

2»7 

2xya' — /3 — 2x^/3' — 2xya' — 13 — 2y‘^j3''\ 

„ / da d/3 da d/3 ^ , 

^yi - +^^>2;—-f ?/--2x7, 

dz dz dz dz 


The paramagnetic current density transforms according 
to jp = jp + pVx, 

[r X = [rx (jp-f pVx)]^ 

= [r X jp]^ + [r X pVx]^ 

= [r X jp]^ + [r X p [xAa, yAa, Ay]]^ 

= [r X jp]^ + p [xyAa - xpAa] 

= [r X jp]^. 


-2/3 - (2x2 + 2y2) p'] . 
modification to /3 would 

in the z component of B and B must be unchanged by both Lz and [r x jp]^ are unaffected. 


So, when we are in any gauge of the form of Eq. (B27l, 
Since any modification to /3 would affect the —2/3 term and when we transform between any of these gauges. 
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